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Abstract 

We use Brownian dynamics simulations of a binary mixture of highly charged spherical colloidal 
particles to illustrate many of the implications of the Random First Order Transition (RFOT) theory 
(Phys. Rev. A. 40 1045 (1989)), which is the only theory that provides a unified description of both 
the statics and dynamics of the liquid to glass transition. In accord with the RFOT predictions, we find 
that as the volume fraction of the colloidal particles ^, the natural variable that controls glass formation 
in colloidal systems, approaches (f)A there is an effective ergodic to non-ergodic dynamical transition, 
which is signaled by a dramatic slowing down of diffusion with the translational constant decaying as 
a power law as — >^ 0^. In addition, using the energy metric (a measure to probe ergodicity breaking 
in classical many body systems) we show that the system becomes non-ergodic as 0a is approached. 
The time t*, at which the four-point dynamical susceptibility achieves a maximum, also diverges as 
a power law near (f)A- Remarkably, three independent measures of dynamics(translational diffusion 
coefficients, ergodic diffusion coefficient, as well t*) all signal that at cpA ^ 0.01 ergodicity is broken. 
The translation diffusion constant, the ergodic diffusion constant, and {t*)~^ all vanish as {(j)A — (t>V 
with both (pA and 7 being the roughly the same for all three quantities. Below 0a transport involves 
crossing suitable free energy barriers. In this regime, the density-density correlation function decays 
as a stretched exponential {exp{—^Y) with (3 0.6. The 0-dependence of the relaxation time is 
well fit using the Vogel-Fulcher-Tamman law with the ideal glass transition occurring at (pK ~ 0.47. 
By using an approximate measure of the local entropy (53) we show that below 0a the law of large 
numbers, which states that the distribution of S3 for a large subsample should be identical to the 
whole sample, is violated. The variations in statistical properties of observables between subsample to 
subsample results in dynamical heterogeneity, which we show is a consequence of violation of law of 
large numbers. The comprehensive analyses provided here for Wigner glass forming charged colloidal 
suspensions fully validate the concepts of RFOT theory. 



2 



I. INTRODUCTION 



The liquid to glass transition is subject of great interest because concepts developed in the 
studies of the structural glass transition (SGT) manifest themselves in other areas in condensed 
matter physics and biology. A growing number of experimental, theoretical, and simulation 
studies have established that the Random First-Order Transition (RFOT) |l| is the most viable 



theory of the SGT currently available. An early review can be found in [2|, and a number 
of recent reviews have provided further analyses, limitations, and applications of RFOT [3-7]. 
The RFOT theory, which explains most of the universal aspects of glass forming materials, 

utions of the statics and dynamics of certain infinite ranged 



was initially inspired by exact so 

or mean-field spin glass models IsHJJI. In the spin glass (SG) models, randomness is explicitly 
introduced in the Hamiltonian by assuming that the interactions between the spins are randomly 
distributed. 



In a subsequent study, Kirkpatrick and Thirumalai [13| (KT) explicitly showed that the 
same general scenario also emerges using equilibrium and dynamical solutions of regular density 
functional Hamiltonian (DFH), appropriate for liquids, thus explaining how randomness is self- 
generated as a liquid undergoes SGT. Here, KT also introduced the idea or technique of inserting 
an infinitesimal replica symmetry breaking field into non-random Hamiltonians so that the 
multiplicity of states in these systems could be investigated. The study in [3] was the first 
paper to show that the ideas that were developed to make precise connections between statics 
and dynamics in SG models, which appeared to be restricted to systems with quenched disorder, 
could be used in any system where there are many statistically distributed metastable states. 
The crucial discovery in KT is that below the dynamical transition temperature, Ta, there 
are an extensive number, e"^ (A^ is the number of particles and a > 0) of states; the liquid 
becomes trapped in one of the large number of metastable states. This finding allowed us to 
produce a consistent static and dynamical theory of the SGT across the entire temperature range 
using a DFH, which is known to be a reasonable starting point for describing the liquid state. 
Indeed, others have also judiciously exploited the observation that the djTiamical transition in 
these systems can be described using equilibrium theories 0, Q. At a lower temperature, Tk 
(< T4), the number of such states becomes non-extensive, and hence the entropy associated with 
the number of states vanishes signaling the ideal SGT. The RFOT theory has been subsequently 



elaborated and expanded on by numerous authors (see for example js, 5, 6.ll6l-ll8| and references 
therein) in an effort to explain glassy states of model systems. It is worth emphasizing en passant 
that in more recent discussions of RFOT only activated processes leading to transport below T4 
are emphasized. Such discussions have created an erroneous impression that RFOT is largely a 
thermodynamic theory whereas the essence of RFOT is to cover the equilibrium and dynamical 
aspects above and below Ta in a consistent manner. Finally, it is worth noting that RFOT is 
qualitatively different from the predictions where the SGT is a regular first order phase transition 



19|. 



Within a mean field picture the barriers separating the multiplicity of disjoint states that 
exist between Tk and Ta increases with the system size so that at the T = T4 there is a genuine 
ergodic to non-ergodic transition. In systems containing particles interacting via short range 
interactions we still expect that for T < Ta there will be finite domains with properties similar 
to the global metastable states predicted by mean-field theory, except that the these states are 
no longer truly disjoint. In this case, ergodicity is effectively broken because the observation 
time scale, Tobs is less than the relaxation time scales, a notion that can only be understood 
using a consistent static and dynamical approach. The long time dynamics below Ta would 
then be governed by activated transport, as the domains change from one metastable state to 
another. Within RFOT the driving force for transport is purely entropic as the system can 



access a large number of states by making transitions between the so-called mosaic states 
The exponentially large relaxation times in glasses are most naturally obtained by considering 
activated transitions between different mosaic states. Very general arguments [li] lead to the 
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conclusion that these droplets are characterized by a (diverging) length scale ^ ~ r~d and a 
characteristic free energy barrier AF'^ ~ where r = ^^^^ , separating two adjacent mosaic 

2 

states. In three dimensions, in which the present simulations are carried out, ^ ~ r~3 and 
AF* ~ ^2. Finally, because there is a distribution of relaxation times associated with the various 
mosaic states we expect that relaxation of various quantities, such as the density-density time 
correlation function, would exhibit a stretched exponential decay, which we demonstrate here 
for colloidal glasses. 

A number of unique features emerge when a glass is pictured as being partitioned into mosaic 
states as T = T4 whose characteristic temperature-dependent size is denoted by ^i{T). For the 
purposes of discussion we do not consider variations in the sizes of the glassy clusters. (1) Within 



each mosaic state structural rearrangements can be rapid but global relaxation requires activated 
transitions, involving nucleation of new domains. Slow structural fluctuations in glasses occurs 
because of entropic driving forces that enable creation and formation of new glassy clusters 
by fluctuation effects. The effects of fluctuations on a system undergoing random first order 
transition are best characterized as T — ?■ T4 using four-point suceptibility, which was first 



19| to differentiate between regular first 



2l| mostly without attribution 



theoretically discussed by Kirkpatrick and Thirumalai 
order transition and RFOT, and subsequently used by many [20|, 
to the original paper. (2) As the degree of supercooling increases not only does ^i(T) grow but 
also the relaxation time associated with particles that cross domain boundaries exceeds Tots, thus 
leading to broken ergodicity. These physical considerations that are embedded in the droplet 
picture of the RFOT [l| also lead to violation of law of large numbers as the STG occurs . 



In this paper, we illustrate some of the RFOT concepts described above using Brownian 
dynamics simulations of binary mixtures of charged colloidal suspensions, whose phase behav- 
ior is extremely rich and interesting on its own. Although considerable amount of work has 
been done to illustrate the phases of monodisperse colloidal particles much less is known about 
binary mixtures of charged micron-sized particles, which readily form the so called Wigner 
glasses. These fascinating systems form Wigner glasses with finite rigidity, as demonstrated in 
a pioneering study by Lindsay and Chaikin [23j. Subsequent molecular dynamics simulations 
confirmed that these low-density systems form Wigner glasses with ease, and could also form 
BCC-like substitutional crystals 2^ ■ The ease of glass formation in these systems were further 
characterized in terms of localized soft modes to describe the nature of activated transport [24]. 



Similar ideas been used in the context of hard sphere glasses 
renewed interest in the study of charged suspensions 
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25( 1 . More recently, there has been 



27( 1 which apparently manifest (with 



some differences) many aspects of the SGT, revealed in simulations of binary mixtures of soft 



sphere systems |28l-l32[| . mixtures of Lennard Jones particles [28|, l33| with additive diameter and 



non-additive diameter 



34| , and possibly even in hard sphere colloidal suspensions [35| . The ad- 



vantage of studying charged suspensions containing micron-sized particles is that, in principle, 
video microscopy can be used to directly image the trajectories of the particles as a function of 
the volume fraction, (p, thus making it possible to validate some of the predictions made in this 
work. 
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II. METHODS 



suspensions, which form Wigner glasses 



Model. Following our previous work 24 1 , we simulated a binary mixture of charged colloidal 



23|. The systems consist of A^^i colloidal particles with 



radius oi and N2 particles with radius 02. The interaction potential between the colloidal 



particles is modeled by Derjaguin-Landau-Verwey-Overbeek (DLVO) potential, which is 



known 



to describe the interactions between highly charged micron size particles in water |36l-l40| . The 

functional form of the DLVO potential used in our simulations is, 

e^ZiZj f exp [qaj] \ / exp [qaj] \ exp [-qvij] 
1 + qai J V 1 + qaj 



where Z^e is the charge of the macroparticle i, q , the square of the inverse Debye screening 
length is q^ = (Zp + ZkPk), Pk is the number density of the k^'^ species, Z is the valence of 
any added electrolyte, and p is the corresponding number density. Because of effective screening 
DLVO potential the DLVO potential is not as long ranged as the Coulomb potential. Thus, 
we did not find it necessary to use Ewald summation. Instead, we neglected, without loss of 
accuracy, interactions beyond a cut-off distance Tc determined by Vij (rc) = O.OOl/csT. The 
parameters used in the simulations, listed in Table 1, correspond to the experimental system 



23| 



Simulation Details. We performed Brownian dynamics simulations (suitable for colloidal 
suspensions in water) by integrating the following equations of motion. 



V,U{n,---,TN)j^ + pD,Ik(t) (2) 

where n (t) is the position of i^'^ particle, [/ (ri, ■ ■ ■ , f^) is J2ij^j '^ij (i^ij)^ Di is the diffusion coef- 
ficient of the i^^ particle, Ri (t) is the random noise satisfying (^Ri (t) ■ Rj (t')^ = QDi6ij6 {t — t') 
with 5ij being the Kronecker delta and 5{t — t') is the Dirac delta function. Because of computa- 
tional considerations we neglected hydrodynamic interactions, whose effects we intend exploring 
in the future. 

Since the quantities of interest in this study are the two-point density- density correlation 
function and the associated fluctuations (four point susceptibility) the integration step, 5t, 
integration must be smaller than the characteristic time scale, a^/Di, in which the two-point 
density decays. The characteristic distance between particles is = p~^ where p = (^1+^^) 
with V being p-dependent size of the simulation box. The neglect of inertial effects is justifled 
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if 5t is larger than the characteristic decay time, miDi/kBT of the velocity correlation function. 
The values of (A; = 1 or 2) are obtained using Dk = where r] = 0.89mPa ■ s for water. 
With these values the range for 6t turns out to be from 10 ns to 1 ms. We chose 7/is for 6t 
as a compromise between accuracy and computational costs. The use of real times is only for 
estimates, and they may not precisely correspond to experimental times. 



The colloidal system is specified by p and temperature T, which we set to 298K. In contrast 
to the usual simulations of mixtures of soft spheres or Lennard- Jones particles in which glass 
formation is driven by quenching T it is p that controls Wigner glass formation 23] . We use the 
volume fraction, = ^{Nia^ + N2al) as a measure of p. In most cases we measure distances 
in units of as = p^ . The 50:50 mixture used in our simulations corresponds to Ni = N2 = 5000. 



We performed equilibration of the system of 10^ particles by placing them initially at the 
sites of body-centered cubic lattice in a periodic simulation box. The size of box was adjusted to 
obtain the desired 0, and it ranged from (2,500 - 5,000) nm (~ (46 — 92)ai) depending on (p. To 
achieve thermal equilibration, we carried out slow- quenching by controlling the concentration of 
electrolytes following the method used in [ssl. The ratio of the concentration of electrolytes to 
colloids p = p/Z^^ J2k ^kPk was initially set to 5, and the equations of motion were integrated 
for 2 X lO^St in the liquid phase (low 0), and 10^6t in a highly jammed glassy state (high 
0). Subsequently, we reduced p by half, and the simulations were further carried out with the 
reduced p. This procedure was repeated until p reached 5/2^°. After reaching the final value, 
p was set to zero. The protocol used here was needed to relax the system especially at high 
values of 0. We performed the same process was for 0=0.01, 0.02, 0.03, 0.04, 0.05, 0.075, 0.1, 
0.125, 0.15, 0.175 and 0.2. After equilibration data were collected for 10^ time steps. Since this 
time is not long enough to obtain structural relaxation dynamics for > 0.075, we performed 
additional simulations for times ranging from 10^6t to 4 x 10^6t when the system reached high 
density jammed state. For example, at = 0.2, total simulation time was 5 x 10^6t, which 
still proves to be not long enough to accurately extract structural relaxation times at high 
values. We generated 20 trajectories at each and ensemble averages, where appropriate, were 
performed over the trajectories. 
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III. RESULTS AND DISCUSSION 



RADIAL DISTRIBUTION FUNCTIONS: 



The radial distribution functions, and the associated structure factors for various volume 
fractions are show in Fig. la. As has been established in many previous studies for a variety of 
systems the static pair functions, gu{r), guir), and g22{r) do not vary significantly except for 
some very minor changes in the first peak. In addition, there is clear evidence for split second 
peak, which is often found in glass forming systems 4l|]. The longer (compared to systems 
interacting with soft sphere and Lennard- Jones potential) range nature and strength of the 
interaction leads to short range order that persists for several values of even at low volume 
fractions. For the highly charged micron-sized particles studied here the effective hard sphere 
diameter is considerably larger than the bare particle size Thus, although technically 



these can be thought of as low-density low-density systems even the lowest value of would 
correspond to a dense equivalent hard sphere liquid. Apart from this observation, it is worth 
reiterating that the pair functions show practically no change even as is increased by twenty 
fold, hence carry no signature of the SGT. 



DEPENDENCE OF DIFFUSION COEFFICIENTS ON cp: 

The mean-square displacement (MSD) as a function of t, 

1 Nc 

{^rlit)) = —Y.{Ht)-r.{o)r) (3) 
with a = 1 (small sized particles) or 2 (large size particles) are shown in Figs 2A and 2B 

2 

respectively. There are three regimes that are discernible. At short times t < we find 
that (Ar^) increases linearly with t with a slope that is proportional to the short time diffusion 
constant, D^o- As is well known this regime represents caging of a test particle by the neighbors. 
At longer times there is a plateau, whose duration increases as increases. Finally, at much 
longer times the particles undergo diffusive motion, and (Ar^) again increases linearly now with 
Da being determined by interaction with particles. 

The values of Dao, obtained from the slopes of the initial increase in (Ar^(t)) do not change 
significantly as increases (Fig. 2C), because the extent of caging is only weakly dependent on 
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(f). However, the values of Da, calculated from the slopes of (Ar^(t)) at long times (see Figs 2 A 
and 2B), decreases rapidly as increases (Fig. 2D), and nearly vanishes at ~ 0.1, which is the 
first signature of the onset of glassy behavior. The values of for a = 1 and 2 are well fit using 
Da ~ (^ — l)"^-", where in analogy with the expectations based on RFOT we have identified 
Ta ~ The fits, shown in Fig. 2E, yield (pA ~ 0.09 and 'Jd ~ (1.0 — 1.2) depending on the 
particle type a. Three comments about these fits and the plateau duration (Fig. 2E) are worth 
making. (1) The values of 7^) are smaller than what is typically expected based on the mode 



coupling theory predictions. (2) We expect based on RFOT predictions [13|, ll9|], and this is 
demonstrated below from an entirely different perspective, that (pA is the characteristic volume 
fraction at which there is an effective ergodic to non-ergodic transition. (3) The duration of 
the plateau in (Ar^) increases rapidly as exceeds 0.09 (Fig. 2E) further indicating that this 
represents the dynamical transition density. 

^-DEPENDENCE OF RELAXATION OF DENSITY-DENSITY CORRELATION 
FUNCTION AND ACTIVATED TRANSPORT: 

The collective variable that slows down > 0yi is Fg{t), the density- density correlation 
function defined as, 

1 N N 

(^) = 1^ E ^'"^^ E ^"''^"'^'^ (4) 
j=i k=i 

where rj (t) is a position of i*^ particle at time t. The ensemble-averaged isotropic self- 
intermediate scattering function {Fq{t)) is estimated by integrating the ensemble averaged 
{F^{t)), with (■ ■ ■) denoting ensemble average, over space with q = \q\- The plots of the time 
dependence of (Fq-(t)) for = 0.01 and 0.02 in Fig. 3A at g = Qmax = 7^ (where is the location 
of the first maximum in the total pair function calculated using both the particle types) show 
the usual vanishing of {Fq^^^{t)) in the liquid state at long times. In Fig. 3 A the solid lines are 
fits of the simulation data to a stretched exponential function (-Fg,„„^(t)) ~ Cexp{—{^Y) where 
the stretching exponent /3 ~ 0.6 is fairly independent of 0. In the inset to Fig. 3A we show the 
dependence of {Fg^^^(t)) for a range of 0, with the corresponding stretched exponential fits with 
/9 ~ 0.6, which is identical to that found in experiments on hard sphere colloids j43|. When 
exceeds about 0.15 the dynamics is so sluggish that {Fg^^^(t)) does not decay fast enough, 
which is an indication that there is a volume fraction, 0^ ^ (pA at which the system undergoes 



an ideal glass transition. 

The dependence of on (p is shown in the inset of Fig. 3A, with the line being the Vogel- 
Fulcher-Tamman (VTF) fit, 

r D ^ 

T-a ~ TyTFeXp[ ^ _ I. (5) 

By fitting to the VTF equation we obtain tvtf = 0.01s, the fragility index D = 23, (pK = 0.47. 
It is clear that the VTF also provides only a semi-quantitive fit of the entire data set. It appears 
to be a better description of the data at high (p. However, the extremely slow dynamics at 
values of far greater than (pA makes it difficult to obtain numerically converged results for 
Ta- Nevertheless, given that (px ^ 4>a gives us confidence that there is a 0a at which 
the relaxation time associated with the Wigner glass essentially diverges. Because of numerical 
uncertainty associated with computing we cannot rule out the possibility that 6 associated 
with Eq.(l5]) can exceed unity as r = — 1) decreases from about 2, the highest for which 
we can reliably compute r^. 



ERGODICITY BREAKING NEAR 0a: 

As has been emphasized previously the important feature at 0a (Ta in systems where tem- 
perature is the relevant variable) is that there is a break down in ergodicity, implying that time 
averages are no longer equivalent to ensemble averages. Of course, because the mosaic states 
are not truly disjoint ergodicity is only effectively broken, which implies Tabs is much less than 
Ta- In order to determine whether such a transition occurs at ^ 0a we calculated the energy 
metric, which is a general measure for assessing the necessary condition for establishing ergodic 



behavior in classical many body systems 



22 



44| . The energy metric is obtained by first generat- 



ing trajectories using two replicas (different initial conditions) of the system at the same volume 
fraction. For each replica we define the time average value of the energy of the i*'* particle, 

Ei{t)=t-' f Ei{t')At' (6) 

where (s) is energy of particle j at time s, and a labels the replica. The energy metric dap (t) 
is defined as 

dap (t) = E E [E'a; it) - E'^p., it)] (7) 



fc=l i=l 
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where is the number of particles of type k, [t) and E-'^ [t) are the energies of particle j in 
replica a and /3 averaged over time t, respectively. If the system is ergodic on the observation 
time scale (robs) then dai3 (t) vanishes as t — )■ Tots- Thus, when ergodicity is established we expect 
that E'^a;i (tq) = E'^ jj-i (r^) independent of a or /3 or i. This is the situation that pertains to 
the liquid phase. However, if ergodicity is broken, as is expected at ^ 0^, then dap (tq,) — )• C 
iC is a constant) suggesting that the two initial states do not mix on the time scale Tabs- 



As argued elsewhere 



22|, it is the development in time rather than any equal time correlation 



functions that distinguishes a glass from a liquid. It can be shown, using scaling-type arguments, 
that (0) /dajs^t) ^ at long times where the ergodic diffusion constant D]^^ sets the 
approximate time scale in which the two configurations {a and (3) mix. Thus, Ndaf3 (0) /da/3 (t) 
which is extensive in both and Tobs in the liquid phase, remains only extensive in N in the 
glassy phase because Tq, ^ Tobs- 

The energy metric as a function of t is shown in Fig. 4A from which it is cleat that when 
exceeds ~ 0.1, the values of da/s (t) essentially saturate indicating a break down of ergodicity. 
The reciprocal of the energy metric, dai3 {0) /dafiit) (Fig. 4B) clearly increases linearly with 
t at low densities and saturate as increases at long times. From the linear dependence of 
dajsiS^) /dap if) we calculated the dependence of on (Fig. 4C). We find that De ^ ^ 
(nearly vanishes) at ~ 0a, which implies that 0a is the volume fraction at which the time 
and ensemble averages start to deviate from each other [l9|. The dependence of -D^; on can 
be fit using De ~ {<Pa ~ 0)'^^ with 0^ ~ 0.12 and 7^ ~ 1.2. Remarkably, 0^ extracted from the 
dependence oi De coincides with the value of 0^ at which diffusion effectively ceases. Thus, 
0A can be identified with the volume fraction at which ergodicity is broken. 



FOUR POINT DYNAMICAL CORRELATION FUNCTION: 

In order to distinguish between the liquid and glass-like states as approaches (pA it is 
necessary to consider fluctuations in multi-particle correlation functions because there is no 
obvious symmetry breaking as the liquid becomes a glass. Sometime ago two of us [lol suggested 
that fluctuations in a suitably chosen order parameter describe the differences in the nature of 
states probed near (pA (or Ta when a liquid is cooled). The rationale for considering multi- 
particle correlation functions is that the natural order parameter that describes the onset of 
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SGT is the two particle Fq{t), which decays to zero in the hquid phase, and saturates in the 
glassy phase at long times (Fig. 3). Thus, it is the only the fluctuations in Fq{t), which plays 
the role of generalized susceptibility, Xii'llt)^ that can distinguish between the liquid and glassy 
states. Since we first showed X4\Fq{t) (note the notation in [19] is different from the one used 
in subsequent development of this concept) can be used to describe genuine ergodicity breaking 
at Tji in mean field p-spin and Potts glasses a number of studies have used our measure (or 



some variants) to produce evidence for growing dynamical 
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correlation length, and emergence of 



45|-|47|. 



dynamic heterogeneity in a variety of model systems 

The four point correlation function X4|Fg is the ensemble averaged variance in Fg (t), and is 
given by 

XMF,{t) = (FAtf) - {FAt)f (8) 
We calculated X4|f,(^) using a moving time averaging procedure in order to minimize numerical 
errors. We show plots of X'i\Fqit) evaluated at qmax for various values of in Fig. 5. As has been 
observed in a number of previous simulations on different systems the peak in xa, {t) increases 
as (f) increases. The dependence of the time, t*, at which X4|f<, {t*) is a maximum is shown in 
the inset of Fig. 5. Although we are unable to compute X4|-F, (^) accurately for (p > 0.075, the 
changes in t* as changes can be fit to a power law. The details of the fit are in the caption to 
Fig. (5). It is noteworthy that (pA extracted from the fit is essentially the same as in the fit of 
De shown in Fig. 4C, thus establishing that the four-point susceptibility does probe the onset 
of ergodic - non-ergodic transition at (pA- 

Although it is most natural to use fluctuations in Fq{t) to determine the four-point suscep- 
tibility others have considered different less obvious variables. One of these is the total overlap 



function, Q (5, t), defined as [48 1 



N 

n {6, t) = N-'J2co ( |rl (t) - rl (0) | , 6) (9) 

i=l 

where (t) is position of the i*'^ particle, N is number of particles, and u {x, S) is step function 
which is 1 when x < 6. Q{6,t) depends on 6 and it is fixed at O.Sa^. The four-point function 
involving the fluctuations in fl {6, t) is defined as, 

x,\n{t) = N[{n{tf)-{n{t)f]. (10) 

In Fig. 6 we show the time evolution of X4|q if) for various (p values. Both Xi\^ if) and X^Fqf) 
evaluated at gmax are nearly identical. The dependence of t* on (p calculated using X4-\q. f) also 

12 



shows a power law dependence (Fig. 6). The only difference is that the exponent characterizing 
the divergence of t* as — ?■ 0a is 1.20 as opposed to 1.05 obtained in Fig. 5. It is gratifying 
that the two independent ways of computing four-point functions yield (nearly) quantitatively 
identical results. More significantly, the values of (f)A as well as the exponents (7/), 7^;, and 
7jx) characterizing the dependence of the translational diffusion coefficients, ergodic diffusion 
coefficient, and t* are similar. 



DYNAMICAL HETEROGENEITY IS A CONSEQUENCE OF VIOLATION OF LAW 
OF LARGE NUMBERS: 

A tenet of statistical mechanics is that many body systems obey the law of large numbers, 
implying that equilibrium properties of a large subsample is on an average identical to the entire 
sample. In the liquid phase (0 < 0a) the statistical properties of any subsample should coincide 
with that of the entire sample provided the subsample contains a large number of particles 
observation time is long enough compared to -D^^. In contrast, in the glassy phase we expect 
that statistical properties (distribution of energies of individual particles for example) of even a 



large subsample can deviate from that of the entire sample 22[. Moreover, one would expect that 
two distinct subsamples, which become equivalent in a liquid when viewed over a short period 
time, would remain inequivalent (or do not exchange) at volume fractions less than 0a even when 
T'ofes-D^^ ^ 1. Thus, no single subsample can characterize the distribution of observables of the 
entire sample in highly supercooled liquids. In other words, in the glassy phase the law of large 
numbers is violated, and there are subsample to subsample fluctuations. Only by examining 
the entire sample on time scale Tobs ^ can the equivalence of time and ensemble averages be 
established. These arguments suggest that the so called dynamical heterogeneity, which is one 



of the characteristic of glass forming systems 49|, |50|] is a consequence of the emergence of glassy 



clusters, which remain inequivalent even on t ^ Tobs- Because of the likely variations in both 
equilibrium and relaxation properties from subsample to subsample a glassy phase is inherently 
heterogeneous in a dynamical sense. 

In order to illustrate the violation of large numbers we first consider an approximate measure 



of structural entropy S3 given by |51| . 



S3 = I y 4vrr2 [g (r) In^ (r) - {g (r) - 1}] dr (11) 
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where g (r) is pair-correlation function, and p is the number density. Parenthetically ye note 



that the entropy density involves the full hierarchy of molecular distribution functions 
we define a local structural entropy measure 53''' for particle j using, 



Here, 



4^) = I / 47rr2 



gii) (^) In ^0) (r) - {(/(j) (r) - l}] (12) 



where ^'3"'^ (r) is pair-correlation with respect to the j*'^ particle. The time evolution of the 
distribution of P{s^\tA) of the time-averaged S3 (see below) for the entire sample and a subsample 
of size ^ ~ 3.3as (containing large enough particles) are used to demonstrate the violation of 
large numbers. The time-averaged local structural entropy associated with particle j is given 

by, 

sHtA)=t-/ r si{s)ds. (13) 







In Fig. 7A we show the distribution, P{s'i\tA) at different values of Ia for = 0.02. As 
increases P{s'i\tA) converges and the its width narrows as expected for a system approaching 
equilibrium. In contrast, P{s^\tA) for cf) = 0.2 (Fig. 7B) is essentially frozen in time indi- 
cating that the transport of particles required for ergodicity to be reached does not occur on 
tA = 12.5De\ It is more instructive to simultaneously compare the time evolutions of a large 
subsample and the whole sample for (p = 0.02 and (p = 0.2. The four panels in Fig. 8 vividly 
show that for the liquid case (0 = 0.02) the distributions P(s3|t^) are almost the same for all 
tA values as is to be expected based on the law of large numbers. In contrast, at higher vol- 
ume fractions (> 0^) where ergodicity is effectively broken, the P{s3\tA) for the subsample are 
substantially different from that of the entire sample, thus violating the law of large numbers. 
Interestingly, there are subsample to subsample variations in P(s3|t^) even with tA = 12. 5D^^ 
as shown in the inset in Fig. 8D. Because different subsamples behave in a distinct manner 
and do not become equivalent the dynamics below is (pA is heterogeneous. Thus, emergence of 
dynamical heterogeneity is due to the violation of large numbers. 

The inequivalence of particles in distinct regions has been probed by assessing how a subsam- 
ple dynamically approaches the full sample using mult i- dimensional NMR |52i]. By measuring 
correlated rotational jump motions it was argued that time for exchanging between different 
subsamples indeed increases dramatically. While not conclusive, these experiments support the 
arguments relating violation of large numbers to dynamical heterogeneity. 

Pictorially we can see how the frozen dynamics is manifested in Wigner glasses. In the top 
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panel in Fig. 9 we show the time evolution of particles within ^ ^ 3.3as for (j) = 0.02. It is clear 
that there are dramatic changes in the configuration, which explains how ergodicity is established 
by particles of a given type becoming equivalent on t ~ D^^. It is worth noting that movement 
of these particles at = 0.02, which is far above 0^, is highly correlated, and appear to undergo 
string-like motion. In sharp contract, the particles at high density (0 = 0.2) are frozen. These 
represent low entropic droplets with considerable local order, which do not propagate across 



the entire sample. The emergence of 
studies on entirely different systems 
(below) (pA (Ta). 
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ocal order, which is consistent with a number of recent 



53| is inconsistent with string excitations at (T) above 



IV. CONCLUSIONS: 

In this paper, we have used simulations of a binary mixture of charged colloidal liquid, which 
forms a Wigner glass, to confirm and relate many aspects of the RFOT of the STG transition. 
In most studies the SGT is studied by supercooling a liquid by varying temperature rather than 
controlling the density, which is the natural variable in colloidal suspensions. In particular, the 
simulations show that dynamics slows down before (pA is approached. Remarkably, three ways 
of measuring the time scale associated with dynamics (translational diffusion coefficient, ergodic 
diffusion constant, and the dramatic increase in the time at which the four point susceptibility 
has a maximum) all exhibit singular behavior ai (p ^ cpA- These results show that cpAi which 
should be identified with the mode coupling cpA-, does represent the onset of breakdown of 
ergodicity, as predicted by considering by the dynamical aspects of p-spin and Potts glasses 
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191]. Although our goal is not to obtain precise values of the exponents 'jD,'jE,a-nd'j^ 
it is gratifying to note that they are approximately equal, which further underscores that the 
avoided dynamical transition density is physically relevant for the transition to the glassy state. 
It is worth stressing that MOT, which can be approximately related to equilibrium static density 
functional theory 5J], does not provide a description of transport above (below) (pA (Ta). Below 
(pA, the time scale for density relaxation increases dramatically, which is intimately connected 



to a growing correlation length 



55 



571]. In particular, the dependence of relaxation times is 
more consistent with the VTF law. The large value of r ^ (^ — 1) prevents us from extracting 
the characteristic exponents that enter the description of transport based on the droplet picture 
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below (j)A PJ- The value of ~ S.Sa^, with r ^ 1.4, is remarkably similar to the estimate 



made sometime ago for mixtures of Lennard- Jones particles 



331]. The consistency of the data 



with VTF does suggest that the length scale, ^, associated with the mosaic states must grow 



35, 



48|, 



57|, 



58j as the (or T) increases (or decreases) diverging near (px or with an exponent 
J/ = I in three dimensions. The value of = | implies that the characteristic barrier height 
separating two mosaic states must scale as AF* ~ ^2 Several simulations have confirmed 
our predictions 53, 59|. Our work further shows that the RFOT concepts [l| linking statics and 



dynamics [13|, which have remained virtually unchanged since its introduction, is also valid for 
charged colloidal suspensions. 

We conclude by briefly discussing the important features of RFOT, and comparing them to 
Adam-Gibbs theory of the glass transition that has a superficial resemblance to RFOT. The 
hallmarks of the RFOT theory of the structural glass transition are: 

1. There is a temperature or density, T4 or ipA-, where there is an avoided mean-field phase 
transition. Around this temperature there is a dramatic, power law singularity slowing 
down in the dynamics of the liquid. For T < T4, or > 0^, the transport involves 
transitions across free energy barriers, i.e., it is activated, hence the subscript A. The 
avoided transition temperature or is identified with the mode-coupling glass transition 
temperature or 0^. 

2. For T < T4 the driving force for the activated transport is a local complexity, which is 
similar to a configurational entropy, Sc-, that vanishes at a lower temperature denoted by 
Tk . For T ^ the average relaxation time in the liquid is exponentially slow and is 
shown to be given by the VFT law. For T < Tk large scale transport ceases to occur. 

3. As T — i- Tk there is a length scale associated with the characteristic size of the glassy 
droplets, ^, which are constantly transitioning from one metastable glassy state to another. 
This length diverges as ^ ~ (T - Tk)"^/'^ according to very general arguments [l|. All 
three features are needed to understand numerous experimental and simulation results. 

All three features, which require a unified treatment of the static and dynamical aspects of the 
SGT, are needed to understand numerous experimental and simulation results. 

The RFOT theory of the SGT can be compared and contrasted to the Adams Gibbs theory. 
Most significantly, in the AG theory there is no such thing as a distinct T4. The transport is 
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always activated, vanishing at a Tr- where some sort of configurational entropy vanishes. The 
slow relaxation time is given by a VFT law. Therefore, it does not include the phenomena 
(ergodicity breaking and emergence of dynamic heterogeneity) associated with T4 or 0^, that 
are illustrated here using Wigner glasses as examples, and are generally associated with the 
MCT of the glass transition. Indeed, the absence of Ta or (pA in any theory of glasses should be 
considered a major weakness because numerous experiments clearly show that the very nature 
of transport changes at these temperatures or densities, perhaps in an universal manner |60 |. 



Second, the AG characteristic length scale behaves as ^ ~ (T — Tk) ^^'^ as T — )■ Tk- The ^ (| 



53 



59|. 



in three dimensions) exponent is not in agreement with simulation data 
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TABLE I. Values of the number of particles (Nk), electrostatic charge in units of e (Zk), radius (afc), 
mass (mfc) and diffusion coefficient (Dk) used in the simulations are presented. Here k refers to particle 
type. The values of are computed using the Stokes-Einstein formula. All parameters are chosen to 
simulate the experimental system of Lindsay and Chaikin [23 1. 





Nk 


Zk 


ak 


nik 


Dk 


k=l 


5000 


300 


545A 


4.3 X 10^ amu 


4.53^m^/s 


k=2 


5000 


600 


iiooA 


34.4 X 10^ amu 


2.24^mVs 
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Figure Captions 

Figure 1: Pair correlation functions for binary charged colloidal particles. gii{r) (A) Radial 
distribution functions, g'ii(r), between type-1 (smaller-sized) particles. The corresponding struc- 
ture factors Sii{q) are in the inset. (B) Same as (A) except these curves are for the larger-sized 
particles. The inset contains 5*22 (?) plots, between type-2 particles, and 5^12 (r) (C) Pair func- 
tions between type-1 and type-2 particles and the corresponding Si2{q)- In all panels distance 
is measured in units of and the scattering wave vector is in unit of Volume fraction, 0, 
presented here are 0.2, 0.175, 0.15, 0.125, 0.1, 0.075, 0.06, 0.05, 0.04, 0.03, 0.02 and 0.01 from 
top to bottom. 

Figure 2: Mean-squared displacements (MSD) as a function of time, t. (A) Time-dependent 
changes in MSD, (Ar^(t)) as a function of type-1 particles. (B) Same as (A) except the plots 
correspond to (Ar|(t)). In both (A) and (B) the values are 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 
0.075, 0.1, 0.125, 0.15, 0.175 and 0.2 from top to bottom. We have normalized, {Al{t)) by a^. 
(C) Diffusion coefficient in the short time limit for type-1 particles (red squares) and type-2 
particles (blue circles) as a funcion of 0. (D) The values of the diffusion coefficients obtained 
from the long time values of (A^(t)) as a function of (f). (E) Fits of diffusion coefficients in (D) 
to power-law ^ (^(f)~^ — 0^^)''^°. The fits yield (f)A — 0.06 for type-1 particles (red, filled square) 
and 0.064 for type-2 particles (blue, filled circle). The values of ^/d are 1-0 and 1.2 for type-1 
and type-2 particles, respectively. (F) Length of plateu In (r^) in (A) and (B) are plotted as a 



function 0. Plateau length is defined by the region where 



din Ar^ 



is less than 0.3. 



dint 

Figure 3: Characteristics of the scattering function. (A) Self-intermediate scattering func- 
tion Fq (t) a,t q — q^nax are presented for various 0. To illustrate the goodness of the stretched 
exponential fits we show Fq^^ (t) at 0=0.02 (blue) and 0.01 (red). The relaxation times Tq, as a 



function of are shown in the inset. The line, a fit to the VFT-law tq exp 



D 



with the 



.(</'a/<^-1) 

exclusion of = 0.175 and 0.2, gives tq = O.Olsec, D = 22.7 and (pK = 0.465. (B) The set 
of curves for the time-dependence of Fg (t) at g = ^max corresponding to = 0.2, 0.175, 0.15, 
0.125, 0.1, 0.075, 0.06, 0.05, 0.04, 0.03, 0.02, 0.01 from top to bottom are displayed. All the 
displayed Fq (t) at q — gmax curves are fit to a stretched exponential function C exp — (t/Ta)^ 
with /3 — 0.6. The quality of the fits justifies the assumption that ^ is a constant. 

Figure 4: Probes of ergodicity breaking. (A) Energy metric, dap{t) as a function of t for 
various 0. The values of 0s presented are 0.2, 0.175, 0.15, 0.125, 0.1, 0.075, 0.06, 0.05, 0.04, 
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0.03, 0.02 and 0.01 from top to bottom. (B) The reciprocal d{Q)/d{t) plots are fit to DeI + Ce 
to extract De (C) Dependence of De on is plotted. The dependence of De on (f) can be fit 
using De ^ {4>a — 4>V'^ with 0.10 and 7^ ^ 1.2. (D) The same data in log-log scale is used 
to illustrate the dramatic decrease in as approaches (j)A- 

Figure 5: Four-point susceptibility Xiif) function determined by fluctuations of self- 
intermediate scattering function Fq (t) for all pairs of particles. We fix q is fixed at the first 
peak of total structure factor calculated irrespective of particle identity. We present results for 
0.02, 0.03, 0.04, 0.05, 0.06 and 0.075 from left to right. The convergence at higher values of 
(f) arc too slow to be reliable. In the inset the position of the peak t* as a function is shown as 
solid squares. The solid line represents a power-law (^(1)^^ — (p^^^ "'^ the data. The fit gives 
(pA ~ 0, 1, and is the is MCT volume fraction for the Wigner glass. Note that this is the same 
value extracted by fitting the dependence of the diffusion coefficients and the ergodic diffusion 
coefficient on 0. The value of the exponent 7^ is 1.05. 

Figure 6: Evolution of four-point function X4|n(^) defined by variation of total overlap 
function Q {S,t) for all pairs of particles. We show x^init) for (f) 0.02, 0.03, 0.04, 0.05, 0.06 and 
0.75 from left to right. Inset shows the dependence of the position of the peak t* in X4\n{t) 
(squares). The solid fine is a fit to a power-law — 0^^) . The fit gives 4>a ~ 0.1. The 
value 7^ is 1.20. 

Figure 7 : (A) Distribution of time-averaged S3 for liquid (0 = 0.02) at various time intervals 
are plotted. Here, time interval is 12.5, 7.5, 5, 2.5, 1.25, 1 and 0.5D^^ from top to bottom. 
(B) The same graph for glass (0 = 0.02) with the same time interval as in (A). 

Figure 8 : Dependence of P{s3\tA) for different values of tA- (A) The red and green curves 
correspond to the entire sample and a subsample with size ^ ~ 3.30^, respectively. The volume 
fraction is = 0.02. The blue curve gives P{s3\tA) for a subsample of a glassy state, and the 
black is the corresponding result for the entire simulation box. The value of tA = 0.5De^. (B- 
D) same as (A) except the values of tA change as indicated. In the inset in (D) we also show 
-^(ssI^a) for another subsample in red. The structural features of the two subsamples are shown 
in Fig. 9. 

Figure 9 : Illustration of time evolution of particles within ^ = 3.30^. The panels on top 
are for a hquid (0 = 0.02) at two times {t = and t = 2De^). Small particles are colored in 
yellow and large particles are in violet. An unmistakable correlated string-like motion is evident 
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especially for the motion of the large particles. The bottom panels show how two different 
subsamples at = 0.2, whose P^s^I (tA)) for = 12. 5D^^ are shown in the inset in Fig. 8D, 
evolve over time. Blue spheres represent small particles and red corresponds to large particles. 
While snapshot shows large difference between (left, top) liquid at initial time and (right top) 
at t = 2D^^, particles are frozen in glassy state thus there is little change between t = (left, 
bottom) and t = 2D^^ (right, bottom) for glass. Most importantly, there are no string-like 
excitations below 0^. Instead, the snapshots show clear evidence for local order 53 1. 
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